o 
o 

(N 
D 



a 

o 
o 



> 

00 
(N 

ON 

o 
o 



X3 
O 

o 



X 



Diffusion accompanying noise induced transport in frictional ratchets 
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We study the noise induced transport of an overdamped Brownian particle in frictional ratchet 
systems in the presence of external Gaussian white noise fluctuations. The analytical expressions for 
current and diffusion coefficient are derived and the reliability or coherence of transport are discussed 
by means of their ratio. We show that frictional ratchets exhibit larger coherence as compared to 
the flashing and rocking ratchets. 
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I. INTRODUCTION 

The phenomenon of noise induced transport (ther- 
mal ratchets) are being discussed extensively in recent 
years. Ratchets are systems that exploit the nonequi- 
librium fluctuations present in the medium to generate 
a directed flow of Brownian particles in the absence of 
a bias P, Q]. Such ratchet models are found to have 
wide range of applications in biological systems and in 
nanotechnology |3|. 

Several physical models like flashing ratchets Q , rock- 
ing ratchets |a| etc., where potential has been taken to 
be asymmetric in space, have been developed. In these 
models to generate noise induced directional transport 
the nonequilibrium fluctuations need to be correlated 
in time. One can also generate unidirectional currents 
in the presence of symmetric ratchet potentials, how- 
ever it has to be driven by correlated time asymmetric 
force Q. There is yet another class of ratchets, namely 
the frictional ratchets which we consider here, where the 
friction coefficient varies in space 00. In such ratchet 
systems it is possible to get unidirectional currents even 
in a symmetric potential in the absence of a net bias. 
Moreover, external fluctuations need not be correlated 
in time. In the presence of external parametric noise 
the particle on an average absorbs energy from the noise 
source. The particle spends larger time in the region of 
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space where the friction is higher and hence the energy 
absorption from the noise source is higher in these re- 
gions. Thus the particle in the higher friction regions 
feel effectively higher temperatures. Hence the problem 
of particle motion in an inhomogeneous medium in pres- 
ence of an external noise becomes equivalent to the prob- 
lem in a space dependent temperature UM- Such 
systems are known to generate unidirectional currents. 
This follows as a corollary to Landauer's blow torch the- 
orem that the notion of stability changes dramatically 
in the presence of temperature inhomogenieties ■ In 
such cases the notion of local stability, valid in equilib- 
rium systems, does not hold. 

Frictional inhomogenieties are common in superlat- 
tice structures, semiconductors or motion in porous 
media. Frictional inhomogenity changes the dynamics 
of the particle nontrivially as compared to the homo- 
geneous case. This in turn has been shown to give 
rise to many counter intuitive phenomena like noise in- 
duced stability, stochastic resonance, enhancement in ef- 
ficiency etc., in driven non-equilibrium systems [llill^. 

Extensive studies have been done on the nature of 
currents and their reversals, stochastic energetics (ther- 
modynamic efficiency) etc., in different class of ratchet 
models. In contrast, there are very few studies which 
addresses the question of diffusion accompanying trans- 
port in ratchet systems 0, Hit . This is intimately 
related to the question of reliability of transport. Diffu- 
sion infact detriments the quality of transport. In our 
present work we focus on the transport coherence in 
frictional ratchets in the presence of external Gaussian 
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white noise fluctuations. 



overdamped Langevin equation of motion is given by 



Transport of Brownian particles are always associated 
with a diffusive spread. When a particle on an average 
moves a distance L due to its velocity, there will al- 
ways be an accompanying diffusive spread, ff this diffu- 
sive spread is much smaller than the distance travelled, 
then the particle motion is considered to be coherent 
or optimal or reliable. This is in turn quantified by a 
dimensionless quantity, Peclet number Pe, which is the 
ratio of current to the diffusion constant. Higher the 
Peclet number, more coherent is the transport. Peclet 
number greater than 2 implies coherence in the trans- 
port 13]. The Peclet numbers for some of the models 
like flashing and rocking ratchets show low coherence of 
transport with Pe ^ 0.2 and Pe ^ 0.6 respectively. 
Another study on symmetric periodic potentials along 
with spatially modulated white noise showed a coherent 
transport with Peclet number less than 3. In the same 
study a special kind of strongly asymmetric potential 
is found to increase Pe to 20 in some range of physical 
parameters jl4| . Experimental studies on molecular mo- 
tors showed more reliable transport with Peclet number 
ranging from 2 to 6 

In our present work we show that system inhomo- 
geneities help in enhancing the coherence in the trans- 
port depending sensitively on the physical parameters. 
For this we consider a simple spatially symmetric sinu- 
soidal potential. Pe of the order of 3 can be readily 
obtained. The noise strength of the external paramet- 
ric white noise fluctuations play a constructive role in 
enhancing the coherence in transport. As opposed to 
this, temperature (internal fluctuations) degrades the 
coherence in transport. 
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with < f{t) >= 0, and < f{t)f(t') >= 2S{t - 1') where 
< ... > denotes the ensemble average and q the coordi- 
nate of the particle. 

The system is then subjected to an external paramet- 
ric additive white noise fluctuating force ^(i), so that 
the equation of motion becomes 
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with < ^{t) >= and < ^(t)^(t') >= 2rS{t - 1'), where 
r is the strength of the external white noise ^{t). The 
corresponding Fokker-Planck equation is given by |l8l | 
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For periodic functions V{q) and 77(g) with periodicity L, 
one can readily obtain analytical expression for particle 
velocity and is given by [13 
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with the generalized potential ip{q) as 
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and A{q) as 
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II. MODEL: 

We consider the overdamped dynamics of a Brownian 
particle moving in a medium with spatially varying fric- 
tional coefficient 7y(g) at temperature T. Using a micro- 
scopic treatment the Langevin equations for the Brow- 
nian particle in a space dependent frictional medium 
has been obtained earlier 0, 0, . The corresponding 



which inturn determines the effective slope of the gen- 
eralized potential ipiq). Hence the sign of 6 gives the 
direction of current which follows from Eqn^ 

In our present work we have taken the potential 
Viq) = Vbsin(g) and i]{q) = 770 [1 - Xsm{q - 0)], 
< A < 1. The phase lag ^ between V{q) and 77(5) 
brings in the intrinsic asymmetry in the dynamics of 
the system. The effective potential ipiq), 6 and A{q) are 
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obtained from Eqns. El and [7| Following references 
[iH liol , one can obtain exact analytical expressions for 
the diffusion coefficient D and the average velocity v as 



D = 



Iqo L 



I+{x) 



L{1 - exp [-L5]) 



(8) 



(9) 



where /+ (x) and (x) are as given below 
f r 

^\^) Jx-L 



x+L 



1 



I-{x) = exp[-V'(a;)] j exp [V'(y)] (If ) 

L here represents the period of the potential {— 2tt in 
our case). Now, the time taken for a Brownian particle 
to travel a distance L is given asr = L/v and the spread 
of the particle in the same time is given as < (A(j)^ 
2Dt. For a reliable transport we require < (Aq)^ >= 
2Dt < L^. This in turn implies that Pe = Lv/D > 2 
for coherent transport. 



III. RESULTS AND DISCUSSIONS 

The velocity (w), diffusion constant {D) and the Peclet 
number (Pe) are studied as a function of different phys- 
ical parameters. All the physical quantities are taken 
in dimensionless form. In particular, velocity and dif- 
fusion are normalized by (Vo/rjoL) and (Vb/r?o) respec- 
tively. Throughout our work we have set Vq and 770 to 
be unity. Similarly, F and T are scaled with respect to 
Vq and Vq/tiq respectively. 

In Fig. n we have plotted the effective potential tp{q) 
to emphasize this. The effective potential is scaled with 
respect to temperature. The nature of currents (in par- 
ticular the direction of current) can be readily inferred 
from the plot of effective potential ip{q) as a function 
of q for (j) = O.Stt and = I.Stt for various values of 
physical parameters as mentioned in the figure caption. 
For (j) values between to tt the effective potential will 
be tilted downward and hence the current direction is 
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FIG. 1: Effective potential tl^iq) as a function of q for Vo ~ 
f, A = 0.9 and T = 1.0. 



positive. For (j) values between tt and 2tt the current is 
in the negative direction. The barrier heights of the ef- 
fective potential decreases with increase in the strength 
of the external noise F. This in turn will lead to an en- 
hancement in the noise induced velocity which will be 
emphasized in later discussions. 

In Fig. 121 we plot v, D and Pe as a function of the 
phase difference (j) between the periodic functions V{q) 
and 77(g) for a fixed noise strength, F and A, the ampli- 
tude of the periodic modulation of r/(g). We have evalu- 
ated these quantities by numerically integrating Eqns. [51 
& Eqn. El using quadrature methods. All the physical 
parameters are mentioned in the figure caption. As ex- 
pected, all the quantities are periodic functions of phase 
6 with the velocity v being zero at = , tt and 27r 
Hi IE 1^ • In the range between to tt the effective 
potential ■)/'(<?) will be tilted down in the forward direc- 
tion and the current hence is positive. When the phase 
difference is between tt and 2 7r the effective potential 
will be tilted in the opposite direction and hence the 
current is in the negative direction as discussed above. 
This also follows from the fact that in the region of (f> 
between to tt higher friction regions lie between the 
potential minima and the corresponding neighbouring 
maxima on the right and hence the particle on the aver- 
age gains more energy in this region as discussed in the 
introduction. As a consequence the particle in a well 
finds it easier to cross the peak of the potential and go 
over to the right side of the well as compared to crossing 
over to the other side (left of the well). Hence current 
in the positive direction is assured. It should also be 
noted that when current is zero diffusion constant is fi- 
nite (at the values of </> = 0,7r and 27r) 0, H 0- For 
the parameters chosen it appears as if v and D are sym- 
metric around 7r/2. However it is not so. This has been 
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verified with other parametric values. Also, between 
and TT and vr and 27r peclet number exhibits a maxi- 
mum with the value being around 2.5. Both the velocity 
and diffusion constant exhibit a maxima between and 
TT. However, the range of variation of the amplitude of 
velocity is higher than that of diffusion resulting in a 
maxima in peclet number. 
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FIG. 2: Plot V, D and Pe vs 
T = 0.01 



for r = 1 , A 



0.9 and 



To understand the role of spatial asymmetry in the 
potential we have considered an asymmetric potential of 
the form V{q) = Vo[sm{q) — sin(2g)] with the asym- 
metry parameter fi ~ 1. Fig. shows the behaviour of 
V, D and Pe for this simple asymmetric case with the 
other physical parameters kept the same as in Fig. [5] 
For this case we notice that the velocity is not zero 
when (j) — 0,TT or 27r as anticipated. Moreover, it is 
also clear that the presence of this simple asymmetry 
does not make a significant contribution to increase an 
increase in Pe. Hence we restrict to the simple poten- 
tial V = Vq sin(a:) in our further analysis. From these 
two figures we conclude that for a certain range of (j) 
depending on system parameters the transport is coher- 



ent (Pe ^ 2) and moreover this range is quite large. In 
Fig. 21 we plot V, D and Pe as a function of the exter- 
nal white noise strength F. The parameters chosen are 
given in the figure captions. The value of (j) corresponds 
to the value at which a maximum in Pe is seen for the 
parameters chosen in Fig. El We observe that the ve- 
locity increases monotonically with the external white 
noise and saturates at higher values of F (not shown 
in figure). In contrast, the diffusion coefficient keeps 
on increasing indefinitely. The Pe exhibits a maxima 
around F = 1.8 and the value being approximately 2.7. 
The reason behind this maxima is same as that given in 
jl3 | for the symmetric potential case. 

In Fig.|Slwe plot velocity, diffusion and Pe as a func- 
tion of r for F = 1.7. We observe that the current 
exhibits a peak and monotously decreases to zero as 
a function of T. At higher T, (fc^T > F770) this be- 
haviour is expected as the temperature overshadows the 
effect of inhomogeneity in the effective potential thus 
suppressing the ratchet effect. If F is much larger than 
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FIG. 3: Plot of V, D, and Pe vs </> for /i = 1, F = 1, A = 0.9 
and T = 0.01 
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FIG. 4: Plot of V, D and Pe vs F for = 0.48 tt, A = 0.9 
and T = 0.01 



FIG. 6: Plot of V, D and Pe vs A for 
r = 0.01 



: 0.48 TT, F = 1.7 and 
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FIG. 5: Plot V, D and Pe vs T for = 0.48 tt, F = 1.7 and 
A = 0.9 



1 the peak in the current disappears and hence cur- 
rent monotonously decreases with T . In contrast, diffu- 
sion constant increases monotonously. The Pe decreases 
monotonously with T thereby suppressing the coherence 
in the transport. 

In Fig. 1^ we plot velocity, diffusion and Pe for values 
oi (f) ~ 0.487r and F — 1.7 as a function of the amplitude 
of oscillation of the friction coefficient A (0 < A < 1). 
We see that the current increases as a function of A 
whereas the diffusion constant decreases as a function 
of A. In the present case diffusion constant exhibits a 
minima near A = 0.8. However, for other parameter 
values we see that D need not exhibit a minima. It 
monotonically decreases. It should be noted that A can- 
not be greater than or equal to 1 so as to maintain the 
friction coefficient to be positive. The peclet number 
increases monotonously as a function of A. This is the 
only case for which we observe a decrease in D while cur- 
rent simultaneously increases as a function of A. Thus 
from our plot we see that the presence of A makes the 
transport more coherent. 
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IV. CONCLUSIONS 

We have studied the rehabihty or coherence of trans- 
port in frictional ratchet systems in the presence of ex- 
ternal white noise fluctuations. Peclet number greater 
than 2 can be readily obtained in a wide range of pa- 
rameter space implying that the transport is coherent. 
External noise (T) helps in transport coherence whereas 
the internal noise (T) degrades the transport coherence. 
As a function of F and T higher drift velocity is also 
linked with higher diffusion coefScient. It is always de- 



sirable to have larger accompanying current with mini- 
mal diffusive spread to have reliable transport. For this 
we should have regions where velocity increases accom- 
panied by a decrease in diffusion coefficient as a function 
of F andT IllH- 

This maybe achieved in our model 
by appropriately choosing strongly asymmetric spatially 
periodic potential like the one used in Q|. Our fric- 
tional ratchet systems exhibits much larger coherence 
compared to previously studied single particle flashing 
and rocking ratchet systems. 
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